First-Passage Exponents of Multiple Random Walks 
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We investigate first-passage statistics of an ensemble of N noninteracting random walks on a line. 
Starting from a configuration in which all particles are located in the positive half-line, we study 
S n (t), the probability that the nth rightmost particle remains in the positive half- line up to time 
t. This quantity decays algebraically, S n (t) ~ t~^ n , in the long-time limit. Interestingly, there is a 
family of nontrivial first-passage exponents, /3i < P% < • • • < Pn—i\ the only exception is the two- 
particle case where /9i = 1/3. In the N — > oo limit, however, the exponents attain a scaling form, 
Pn{N) — > P(z) with z = (n — N/2)/y/~N. We also demonstrate that the smallest exponent decays 
exponentially with N. We deduce these results from first-passage kinetics of a random walk in an 
iV-dimensional cone and confirm them using numerical simulations. Additionally, we investigate the 
family of exponents that characterizes leadership statistics of multiple random walks and find that 
in this case, the cone provides an excellent approximation. 

PACS numbers: 02.50.Cw, 05.40.-a, 05.40.Jc, 02.30.Em 



I. INTRODUCTION 



Ensembles of ordinary random walks in one dimension 
are used to model physical, chemical, and biologi- 
cal processes ranging from wetting p| M to the motion of 
colloidal particles in narrow channels [7 , H| and reaction- 
diffusion processes @- In particular, first-passage prop- 
erties (To| of multiple random walks explain dynamics of 
interacting spins [TJ Ell > l ea d changes in a voting process 
0, [H[ , and lifetime of knots in polymer chains [14| . 

In these examples, first-passage properties are inter- 
twined with the ordering of the walkers [|| [l5l - |22j |. 
Previous studies focused on first-passage statistics of ex- 
treme particles; for example, the probability that a single 
prey particle survives the predators to its left [l|| Il8j . 
These problems involve a single first-passa ge p robability 
and hence, a single first-passage exponent [1(| . 

In this paper, we ask first-passage questions that con- 
cern the bulk particles, not necessarily the extreme ones. 
We find that a family of first-passage exponents char- 
acterizes the first-passage kinetics. These exponents de- 
pend on two parameters: the particle order and the total 
number of particles. Yet, when the number of particles 
is very large, the exponents depend on a single scaling 
variable. This scaling behavior is unusual. In equilib- 
rium as well as in non-equilibrium settings, one or two 
scaling exponents quantify a scaling behavior [H, [24j . 
In the present case, however, the exponents themselves 
obey scaling, and remarkably, there are scaling laws for 
the scaling exponents. 

We consider an array of N identical particles that un- 
dergo simple random walk on a one-dimensional line, and 
investigate two different first-passage problems where the 
particle order plays an essential role. The first problem, 
discussed in Sections II-V, concerns the probability that 
the nth rightmost particle hasn't crossed the origin, if all 
particles start to the right of the origin (Figured]). This 
quantity equals the likelihood that N — n + 1 out of all 
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FIG. 1: Illustration of a four-particle system. Initially, all 
particles are located in the positive half-line. At time t, the 
position of the rightmost particle is xi(t) and that of the 
leftmost particle is Xi(t). 



random walks have yet to simultaneously reside in the 
negative half-line. In the second problem, described in 
Section VI, we identify the particle order with rank, such 
that the rightmost particle is viewed as the leader, and 
similarly, the leftmost particle as the laggard. We ask: 
what is the probability that the rank of the initial leader 
does not fall below some specified threshold? 

In both problems, we find a family of nontrivial first- 
passage exponents. In both cases, the exponents depend 
on two variables: the particle order and the total num- 
ber of particles, N. Interestingly, in the large-iV limit, 
the exponents become a function of a single scaling vari- 
able. However, the similarities between the two problems 
end here as the two scaling variables are fundamentally 
different and moreover, the two scaling functions are dis- 
similar. 

Our analysis relies on mapping the N noninteracting 
random walks onto a single compound random walk in N 
dimensions. By combining this mapping with exact and 
asymptotic properties for kinetics of first-passage inside 
a cone [2o| . we obtain approximate values for the first- 
passage exponents. The cone approximation is straight- 
forward to implement, yet it yields useful estimates for 
exponents and in particular, this framework faithfully 
captures typical and extremal properties of the first- 
passage exponents. 



II. THE FIRST-PASSAGE PROCESS 

Our system consists of N identical particles. Each par- 
ticle undergoes a random walk on a one-dimensional lat- 
tice. At every time step, one particle is selected at ran- 
dom, and it moves to the left, x — > x — 1, or to the right 
x — > x + 1, with equal probabilities. Time is augmented 
by the inverse number of particles after each such step, 
t — >• t + 1/N. The particles always undergo independent 
random walks and hence, they are noninteracting. 

Let x n (t) be the position of the nth rightmost particle 
at time t (Figured]). We consider the initial configuration 
where all particles are located in the positive half-line 
(Figure [T]), x n (0) > for all n. We stress that it is not 
the initial order, but instead, the order at time t, that 
sets the index n. 

We are interested in S n (t), the probability that the nth 
rightmost particle remains in the positive half-line until 
time t. Hence, S n (t) is the likelihood that x n (r) > for 
all < t < t. In particular, Si(t) is the probability that 
the rightmost particle has yet to cross the origin, while 
Sn (t) is the probability that the leftmost particle has not 
crossed the origin. Clearly, the probabilities S n decrease 
monotonically with n, 



Si(t)>S 2 (t)>--->S N - 1 (t)>S N (t). 



(1) 



The "survival" probabilities S n (t) generalize the classic 
survival probability of a single one-dimensional random 
walk in the presence of a trap. Indeed, when N = 1, we 
have Si(t) ~ i -1 / 2 [10]. As usual, the survival probabil- 
ity S n immediately gives the first-passage probability as 
[— dS n /dt) x dt is the probability that the nth rightmost 
particle crosses the origin for the first time during the 
infinitesimal time interval (t, t + dt). 

The analytically solvable case of two particles yields 
valuable insights into the general behavior. When N = 2, 
we map the two random walks onto a single random walk 
in two dimensions. The position of the two-dimensional 
walk is specified by the positions of the two independent 
walks. At t = 0, the two-dimensional walk is always 
inside the first quadrant (Figure [2]) . 

The quantity Si(t) is the probability that the coordi- 
nates of the two random walks have not become negative 
simultaneously up to time t, or equivalently, the proba- 
bility that the two-dimensional walk remains in the ex- 
terior of the third quadrant (Figure [2^). To find the 
survival probability, we impose an absorbing boundary 
condition along the edge of the third quadrant (Figure 
[5^). Then Si(t) is the probability that the compound 
particle avoids this absorbing boundary up to time t. 
The region in which the random walk can move is a two- 
dimensional cone, or equivalently, a wedge with opening 
angle a — 37r/4. (The opening angle a is defined as the 
angle between the cone axis and the cone surface, so it 
lies within the bounds < a < n.) Thus, Si(t) equals 
the survival probability of a particle diffusing inside a 
wedge with an absorbing surface. This survival probabil- 
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FIG. 2: The compound two-dimensional random walk. The 
line indicates the random walk trajectory and the bullets show 
intermediate locations along this trajectory. Thick lines mark 
the absorbing boundary, (a) For n = 1, the random walk is 
confined to the exterior of the third (negative) quadrant, (b) 
For n = 2, the random walk is confined to the interior of the 
first (positive) quadrant. 



ity decays algebraically, 
S(t) 



t~ 



t/4o 



(2) 



in the long-time limit [loj . By substituting a ~ 37r/4 
into the general expression @, we find the intriguing 
behavior 



Si(i)^- 1/3 , 



(3) 



as t —> oo. As in $Z§, this asymptotic behavior holds 
regardless of the initial position, although the prefactor 
does depend on the initial conditions. 

Along the same lines, S^i) is the probability that the 
positions of the two random walks remain positive up 
to time t, or alternatively, the probability that the two- 
dimensional walk remains in the interior of the first quad- 
rant (Figure (2b)- Now, the boundary of the first quad- 
rant is absorbing (Figure [2Jd), and the random walk is 
confined to a wedge with opening angle a = n/4. Using 
@ , we again find power-law behavior albeit with a larger 
exponent, 



(4) 



As shown below, this particular behavior follows from an 
elementary argument. 



III. A FAMILY OF EXPONENTS 

The analytic results for a two-particle system sug- 
gest that generally, the survival probabilities decay al- 
gebraically, 



s n (t)~t-i>», 



(5) 



in the long-time limit. Moreover, we expect that the 
decay exponents are distinct, and that /?„ = f3 n (N) in- 
creases monotonically with n, 



<(3 2 < ■■■ <Pn-i < Av, 



(G) 




FIG. 3: The survival probabilities Si(t) and for three 

particles. These results are from numerical simulations. 
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TABLE I: The exponent /3 n , obtained from numerical simu- 
lations for N < 6. 

consistent with ([TJ. 

The largest exponent, /3/v, is trivial. The probabil- 
ity 5jv(0 that the leftmost particle has yet to cross 
the origin equals the probability that not a single par- 
ticle crossed the origin. Since the particles are indepen- 
dent, this probability equals the product of the individ- 
ual probabilities for each particle not to cross the origin, 

Sn ~ (t^ 1 ^ 2 ) ~ t~ N ' 2 . Therefore, the largest exponent 
is proportional to the total number of particles, 

Pn = \N. (7) 

Our numerical simulations confirm that indeed, there 
are N distinct first-passage exponents. Figure [3] shows 
the survival probabilities S\ and S2 for a three-particle 
system, and Table I lists the exponents /3 n for TV < 6. 
The simulations confirm that there is a family of N — 1 
nontrivial exponents, {0%, /?2, • ■ • , Ay— i}- 

We performed extensive Monte Carlo simulations us- 
ing the following algorithm. Initially, all N particles are 
located at the same position, x = 1. In each subsequent 
step one particle is selected at random and it moves by 
one lattice site either to the left or to the right with equal 
probabilities. After each step, time is augmented by the 
inverse number of particles, t — > t + 1/N. Throughout 
this random process, we keep track of the total number 
of particles in the region x < 0, and obtain the probabil- 
ity S n (t) from the fraction of realizations in which this 
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FIG. 4: The function t^ 1 Si(t) versus t, using different values 
of Pi. The data is from numerical simulations of a three- 
particle system. 



counter did not exceed n until time t. We computed the 
exponents /3„ by analyzing the local slope — d\nS n / dint 
or alternatively, by seeking the value of f3 n for which the 
quantity t l3n S n (t) achieves a plateau (Figure 2]). 

We verified that the measured exponents f3 n are ro- 
bust. First, we checked that they are not sensitive to 
the initial positions. Second, we used parallel dynam- 
ics where particles move synchronously, rather than se- 
quentially, and found that the exponents do not change. 
The exponents reported in this study represent an av- 
erage over a large number of independent realizations. 
The total number of realizations varies from 10 6 for slow 
first-passage processes (/3 <C 1) to as high as 10 12 for fast 
processes (f3 ^> 1). We note that it is difficult to mea- 
sure tiny exponents f3 <C 1 because the simulations must 
run to very large times, and furthermore, it is difficult to 
obtain large exponents /3 ^> 1, because now, we need an 
enormous number of realizations. 

As a side note, the behavior (JS])-© has a natural in- 
terpretation in the context of single-file diffusion where 
the random walks interact by hard core exclusion [2(l - [32j . 
Under the standard transformation where the identities 
of two particles are exchanged whenever their trajectories 
cross, the noninteracting particle system is equivalent to 
an interacting particle system. In the context of this ex- 
clusion process, the quantity S n (t) is the probability that 
the nth out of N particles avoids the negative half-line 
up to time t. 



IV. CONE APPROXIMATION 

In general, we map the N one-dimensional random 
walks onto a single random walk in N dimensions. To 
find the survival probability S n (t), we require that the 
random walk remains inside the region in which the num- 
ber of negative coordinates is smaller than n up to the t. 
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TABLE II: The exponents /3Jj ono , obtained using the cone ap- 
proximation, that is, Eqs. ((8))- (|ll[) . 



This region occupies a fraction V n of space given by 

N 



Vn = 2- N £ 



N 



(8) 



For example, when ]V=2we have (Vi,V 2 ) = (3/4, 1/4) 
as in Fig.il while (Vi,V 2 , V 3 ) = (7/8, 1/2, 1/8) for N = 3. 

In two dimensions, the random walk is confined to 
two-dimensional cones of opening angles a — 3tt/ '4 and 
a = 7r/4 when n — 1 and n — 2, respectively. In the 
cone approximation, we replace the region of space that 
confines the walk with an unbounded circular cone in N 
dimensions that occupies the same fraction of space V n . 
An A^-dimensional cone with opening angle a occupies 
the fraction of space 



V(a) 



Iode(sm9) 



N-2 



N-2 



(9) 



Using a spherical coordinate system, this expression fol- 
lows from dfl oc (sm9) N ~ 2 d9 where f2 is the solid angle 
and 9 is the polar angle. We now choose the cone opening 
angle a = a(n, N) as follows 



V(a) = V n 



(10) 



The survival probability of a particle that diffuses 
inside an unbounded iV-dimensional cone decays alge- 
braically, S ~ i - ^, in the long-time limit [13, [H, [3J]. 
The exponent /3 = /3(a) decreases as the opening angle 
a increases. In two dimensions, f3 — ir/Aa as in ([5]), and 
in four dimensions (3 = (tt — a) /2a j25[. Generally, the 
decay exponent (3 is the smallest root of a transcendental 
equation involving either P^p+s or Q2P+5, that is, one 
of the two associated Legendre functions 35] of degree 
2/3 + S and order S = ^ 



P. 



(cos a) 



Q2P+& 



N odd, 
N even. 



(11) 



In three dimensions, the transcendental equation involves 
the Legendre function, P2p{cosa) = 0. In all dimen- 
sions, a cone with a = n/2 is simply a plane, and hence 
/3(tt/2) = 1/2. _ 

By construction, the cone approximation is exact in 
two dimensions. For three particles, the cone approxi- 
mation provides good estimates for the smallest and the 



largest exponents. When n = 1, the cone specified by 
([TU]) has opening angle a\ = tt — 0.722734 , and equation 
(fTTj) gives the first-passage exponent /3™ no = 0.216785 
where the superscript indicates outcome of the cone ap- 
proximation. This value represents a very good approx- 
imation, as the simulations yield (3\ = 0.228. For the 
largest exponent (n = 3) the opening angle is a 3 = 
0.722734 and the consequent value ^ ono = 1.407279 is 
close to the exact value f3 3 = 3/2. When n = 2, there is 
a larger discrepancy: 02 — tt/2 leads to the approximate 
value f32° nc = 1/2, while the simulations yield /?2 = 0.68. 

Table II lists the outcome of the cone approximation 
for 2 < N < 6. Comparing with the respective values 
in Tables I, we see that quantitatively, the cone approx- 
imation deteriorates as N grows. Nevertheless, the cone 
approximation faithfully captures all qualitative features 
of the first-passage exponents, as shown below. 



V. SCALING AND EXTREMAL PROPERTIES 

We are especially interested in the behavior when the 
number of particles is large, N 1. Let us first evaluate 
the fraction V n in the limit N — > 00. Using the Stir- 
ling formula N\ ~ \Z2ttN N n e~ N we write the leading 
behavior of the binomial term in ([5]), 



( N ) ~ ^Ne-^ m - N ^l N . 
\mj 



We now take the continuum limit and convert the sum on 
the right hand side of ([8]) into an integral. The quantity 
V n becomes a function of a single scaling variable 



V^JV) -> i erfc 



V2] with z — 



N/2 



N 



(12) 



where erfc(£) = (2/y/n) f^°e~ u2 du is the complemen- 
tary error function. Equation (1121) is valid in the limit 
n, N — > 00 with the scaling variable (n — N/2)/\/N finite. 

Similarly, we evaluate the leading large- N behavior of 
V(a) given in ([9]). The dominant contribution to the 
integral comes from a narrow region of order 1/y/N in 
the vicinity of a = tt/2 where the integrand is Gaussian, 



(sin 9) 



N-2 



-N{tt/2-0Y /2 



Using the leading asymptotic behavior ofthe denomi- 
nator, exp [-N(tt/2 - 9) 2 /2]d9 -> v / 2r7/N, we find 
that the fraction V(a) has the scaling form 



U(a,A0-^erfc(-|j 



with y = (cos a)VN. (13) 



In writing this equation, we used the fact cos a ~ tt/2 — a 
when tt/2 — a — > 0. Equation (fT3| is valid in the 
limit tt/2 — a — > 0, N — > 00, with the scaling variable 
(cosa)v^V finite. Meanwhile, asymptotic analysis of (fTT| 





FIG. 5: The exponent ft versus the scaling variable z. The 
Monte Carlo results are from simulations with TV = 10 3 (bul- 
lets) or TV = 10 4 (squares) particles. The solid line is the 
scaling function (|16[) obtained using the cone approximation. 



shows that in the limit TV — > oo, ir/2 — a ^ with the 
scaling variable (cos a)y/W finite, the exponent P(a) be- 
comes a function of the scaling variable y alone, 

P(a,N) ^ P(y) with y = (cos a)y/N. (14) 

The dependence of the exponent j3 on the scaling vari- 
able y is specified through the transcendental equation 
D2d{y) = 0, where D v is the parabolic cylinder function 
35] of order v, and the acceptable root is the smallest 
one [25j . 

Comparing equation (|12[) with equation (|13j) . we con- 
clude that the first-passage exponent depends on a single 
scaling variable, 



p n (N)->P(z) with 



n-N/2 



TV 



(15) 



in the large- TV limit. Using y = 2z, the exponent p and 
the scaling variable z are related by the transcendental 
equation 



D 2p (2z) = 0. 



(16) 



In particular, P(z = 0) = 1/2 (Figure [5]). 

The scaling behavior (TT5l - p^|) implies that the expo- 
nents p n are of order one only in a window of size y/~N 
centered on the mid-point n = N/2. Otherwise, the ex- 
ponents are exponentially small when N/2 — n ^> y/~N, 
or algebraically large when n — N/ 2 ^> y/N. The limiting 
behaviors of the scaling function captures these extremal 
properties, 



P c ° nc (z) 



\Jz 2 /2it exp(-2z 2 ) z 
z 2 /2 z 



(17) 



Both of these limiting behaviors follow from {To}, see 
[25|. The algebraic behavior in the limit z —± 00 implies 
Pn ~ TV, in agreement with ([7]). 



FIG. 6: The survival probability S(z, t) versus time t. Shown 
are results of simulations with different number of particles, 
TV = 10 4 (bullets) and TV = 10 5 (squares), at three different 

= — 1/2 (upper curves), z — 
1/2 (lower curves). 



values of the scaling variable 
(intermediate curves), and z 



Results of massive numerical simulations with a large 
number of particles confirm (Fig. [5$ that the first-passage 
exponents adhere to the scaling form f|15|). Moreover, the 
shape of the scaling function is qualitatively similar to the 
shape of the scaling function (|16p. 

Additionally, the simulations reveal (Fig. |6|) that not 
only the scaling exponents depend on a single scaling 
variable, but the entire survival probability is also func- 
tion of the scaling variable z. That is, 



S n (t)^S(z,t) with 



n-N/2 



(18) 



when TV 3> 1. We confirmed this behavior numerically 
using a very large number of particles (Figure |BJ). The 
scaling behavior (ITBl is notable because it applies at all 
times: it holds at short times, at intermediate times, and 
at long times. 

Of special interest is the survival probability of the me- 
dian particle which becomes completely independent of 
the total number of particles when TV — > 00. There is an 
interesting limiting value that characterizes the survival 
probability of the median particle, 



?mid = lim Pn/2- 



(19) 



Numerically, /3 m id = 1.16, while the cone approximation 
gives Pmid = 1/2. This universal behavior and the scaling 
form (|15[) show that there is a narrow "scaling-window" 
of width An with An/TV ~ TV~ 4 / 2 centered on the me- 
dian particle. Hence, only the relatively small number of 
particles residing inside this region influence the median 
particle. This behavior is in line the narrow range of cone 
opening angles, Aa ~ TV -1 / 2 , where the first-passage ex- 
ponents change rapidly according to (|14[) . 

The smallest exponent Pi underlying the decay of 
Si(t) is especially intriguing because Si(t) has a sim- 
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TABLE III: The smallest exponent as a function of N. Listed 
are the outcome of the cone approximation, /3l° nc , and the 
results of numerical simulations, ft. 



pie interpretation — it is the probability that the N- 
dimensional walk remains in the exterior of the general- 
ized A^-dimensional "quadrant" (see Figure [5]), a region 
that occupies a fraction V\ = 1 — 2~ N of space. Ta- 
ble III shows that the cone approximation yields use- 
ful estimates for f3i when N is small. The cone ap- 
proximation is also useful for understanding the large- 
N behavior. Substituting V% = 1 - 2~ N into ([10|) and 
[sin(7r - 6)} N ~ (?r - 0) N into ©, we readily find that 
the opening angle approaches a constant, 



(20) 



in the limit N — > oo. In a cone of fixed opening angle 
a > 7r/2, the survival exponent /3 shrinks exponentially 
with dimension, f3 ~ [sin(7r — a)] N in the large- N limit 
[25| . Thus, we conclude that the smallest exponent de- 
cays exponentially with the total number of particles, 

-civ 



ft 



(21) 



with C cono = - ln(sin §) = 0.735167. The simulation re- 
sults are consistent with this exponential decay (Fig. [7]). 

The exponential decay (|21[) and the algebraic growth 
([HI) imply the following limiting behaviors of j3 in the 
scaling regime, 



Pi*) 



exp (—const, x z 2 ) 

z 2 



z — > — oo, 
z oo. 



(22) 



These behaviors are fully consistent with (IT71) . Therefore, 
the cone approximation captures all qualitative features 
of the first-passage exponents including typical behavior 
that includes the form of the scaling variable and the 
shape of the scaling function, as well as extremal behavior 
comprising of exponential decay of the smallest exponent 
and algebraic growth of the largest exponent j3(| . 



VI. LEADERSHIP STATISTICS 

The cone approximation is useful in other contexts. 
We now apply this approach to understand leadership 




FIG. 7: The smallest exponent, ft, versus the total number 
of particle N. The results are from Monte Carlo simulations. 



statistics of multiple random walks. 

We study the very same system of N independent ran- 
dom walks in one dimension. As shown in Figure [TJ 
we rank the particles by position with n = 1 being the 
rightmost particle and n = N being the leftmost parti- 
cle. As if the random walks participate in a competition 
[l8l . [l9l I37rl39| , we view the rightmost particle as leader 
and the leftmost particle as laggard. We now investi- 
gate P n (t), the probability that the rank of the particle 
initially in the lead, does not fall below n up to time 
t. Thus, Pi(t) is the probability that the original leader 
never loses the lead, while Pjv-i(t) is the probability that 
the original leader does not become the laggard. 

If N = 2, we have Pi ~ i -1 / 2 . For three particles, 
Pi(t) - t~ 3 / 4 and P 2 - <~ 3 / 8 19]. In this case, the 
compound random walk in two dimensions is confined to 
wedges formed by two intersecting planes. These wedges 
occupy fractions V\ — 1/3 and V2 = 2/3 of space, and 
the decay exponents are given by 7 = 7r/4a or equiv- 
alent^, 7 = 1/4V. Based on the results for two- and 
three-particle systems, we expect that the probabilities 
P n decay algebraically, 



Pn(t)~t->« 



(23) 



in the long-time limit. Once again, there is a family of 
first-passage exponents, 



71 > 72 > ■ • • > 7w-i- 



(24) 



The largest exponent 71 characterize the probability that 
the initial leader maintains the lead, and the smallest ex- 
ponent 7at_i characterizes the probability that the orig- 
inal leader never turns into the laggard. 

In this case, the compound A^-dimensional walk moves 
in a region that occupies a fraction 



N 



(25) 



of space. For example, to see that V\ — 1/N we note that 
the total A^-dimensional space is divided into N equiv- 
alent regions in which the identity of the leader is the 
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7i 


cone 
72 


cone 
73 


cone 
74 


cone 
75 


3 


3/4 


3/8 








4 


0.888644 


1/2 


0.300754 






5 


0.986694 


0.612148 


0.408397 


0.253371 




6 


1.062297 


0.701508 


1/2 


0.351915 


0.220490 



TABLE IV: The exponents 7™ nc , obtained using a cone in 
TV — 1 dimensions with fraction of space V n = n/TV. 



N 


7i 


72 


73 


74 


75 


2 


1/2 










3 


3/4 


3/8 








4 


0.913 


0.556 


0.306 






5 


1.02 


0.676 


0.454 


0.265 




6 


1.11 


0.769 


0.557 


0.391 


0.234 



TABLE V: The exponent 7„ from numerical simulations. 



same. It is straightforward to generalize this result to all 
n. We have seen that when there are three particles, the 
boundary of the allowed region forms a two-dimensional 
cone. We therefore replace the allowed region with an 
N — 1 dimensional cone occupying the same fraction of 
space (|25|) as the allowed region. The opening angle of 
this cone, a, satisfies 



\JV-3 



N-3 



n 

N' 



(26) 



Given n and TV, we first determine the opening angle a 
by solving (|26p. and then compute the exponent 7 as 
the smallest root of a transcendental equation involving 
one of the two associated Legendre functions with degree 
27 + /i and order fi 



jV-4 
2 



COS OL I 







N even, 
N odd. 



(27) 



By construction, the cone approximation is exact for 
three particles. The cone approximation values, listed in 
Table IV, provide very good estimates, given the Monte 
Carlo simulation values listed in Table V. The numeri- 
cal simulations are implemented using the algorithm de- 
scribed in Section III, except that at time t = 0, the 
random walks occupy TV consecutive lattice sites. 

From (|26|) . we expect that the exponent 7 depends on 
a single scaling variable, 



7„(iV) — > j(x) with x 



n 

N' 



(28) 



in the limit AT — > 00. When TV is very large, we can 
replace the left-hand side in (|26[) with the right-hand side 
in ©. Substituting (|25) into (H3J) shows that 



y 



(29) 



N 


cone 

7jv-i 


7jv-i 


2 


1/2 


1/2 


3 


3/8 


3/8 


4 


0.300754 


0.306 


5 


0.253371 


0.265 


6 


0.220490 


0.234 


7 


0.196216 


0.212 


8 


0.177469 


0.190 


9 


0.162496 


0.178 


10 


0.150221 


0.165 



TABLE VI: The smallest exponent as a function of TV. Listed 
are estimates from the cone approximation, 7^ t \ , and results 
of numerical simulations, 7jv_i. 



N 


_ cone 

7i 


7i 


3 


3/4 


3/4 


4 


0.888644 


0.91 


5 


0.986694 


1.02 


6 


1.062297 


1.11 


7 


1.123652 


1.19 


8 


1.175189 


1.27 


9 


1.219569 


1.33 


10 


1.258510 


1.37 



TABLE VII: The largest exponent as a function of TV . Listed 
are the estimates from the cone approximation, 7i° nc , and 
results of numerical simulations, 71. 



Here erfc _1 (a;) is the inverse complementary error func- 
tion, (see below equation (TT2"j) L Using the scaling behav- 
ior (|16|) . we find that 7 and x are related by the tran- 
scendental equation (Figure [8]) 



D 



2- 



(V2 



erfc" 1 (2x)) = 0, 



(30) 



where D v is the parabolic cylinder function. We note 
that 7(0) = 1/2. Given the definition of x, the scaling 
behavior in the leadership problem is completely different 
than the scaling behavior in the origin-crossing problem. 
Using the asymptotic behaviors 



2-(v / lM 5 )ex P (-a £ 



erfc(0 



we deduce that the variable y defined in (|29|) has the 
asymptotic behaviors 



0, 
1. 




2(1-2;) 



Substituting these expressions into the asymptotic be- 
havior of the scaling function in (|14[) . we find the limiting 
behaviors (see also [25]) 



(31) 




2(1=5) 



8 



y i 




FIG. 8: The exponent 7 versus the scaling variable x. The 
line shows the prediction of the cone approximation given in 
equation (I30p . as well as results of numerical simulations with 
N = 10 3 (bullets) and N = 10 4 particles (squares). 



The exponent 7 decreases monotonically with x: it 
weakly diverges when x — > and it vanishes as x — > 1 . By 
substituting x = l/N and 1 — x = 1/N, respectively, into 
the appropriate expressions in (|31[) . we find the leading 
large- N behaviors 



In 



\\nN 7i = l, 
^lnV n = N-l. 



(32) 



Both expressions match estimates based on heuristic ar- 
guments [18l . ll9| . The smallest and the largest exponents, 
listed respectively in Tables VI- VII, show that the cone 
yields an excellent approximation. Yet, the quality of the 
approximation declines ever so slightly as N increases. 

Monte Carlo simulations with a large number of par- 
ticles confirm the scaling behavior (|28|) . Moreover, we 
numerically verified that the entire survival probability 
becomes a universal function of the scaling variable x, 
in analogy with ([T8"]). that is P n (N,t) — > P{x,t) with 
x = n/N, as N — !• 00. Interestingly, the numerical re- 
sults strongly suggest that the scaling function specified 
in (|3U)1 is asymptotically exact (Figure [5]) . Finally, the 
exponent 7jv/2i that characterizes the probability that 
the original leader always ranks higher than median has 
a simple limiting value, 7^/2 — > h- 



VII. DISCUSSION 

In summary, we investigated two first-passage prob- 
lems involving ordered random walks in one dimension. 
In both cases, there is a series of survival probabilities 
and a nontrivial family of first-passage exponents. A uni- 
versal function describes the exponents when the number 
of particles is very large. Remarkably, there are scaling 
laws for the scaling exponents. 

In general, a first-passage process with N random 
walks is equivalent to diffusion in a high-dimensional 



space with a complicated absorbing boundary. This 
boundary is typically formed by multiple intersecting 
planes. Solving the first-passage problem is equivalent to 
solving the electrostatic problem given these boundary 
conditions [25l. IZcil - II^ . For example, to find /3i = 0.228 
for a three-particle system, we must solve for the poten- 
tial in the exterior of an insulating three-dimensional cor- 
ner. Yet, the solution, which likely requires an ingenious 
implementation of the image method, remains unknown 
even in this physically relevant geometry [43| . 

To circumvent this difficulty, wc introduced the cone 
approximation in which the absorbing boundary is re- 
placed with the surface of a suitably chosen cone. This 
approximation is straightforward to implement and uti- 
lizes exact and asymptotic properties for first passage in a 
cone. The cone approximation provides useful estimates 
for the exponents and moreover, it faithfully captures 
both typical and extremal features. In particular, the 
cone approximation yields the correct scaling variable for 
the first-passage exponents. 

We used the term "cone approximation" , yet this 
framework produces the lower bound whenever the al- 
lowed space is a generalized cone with the crucial prop- 
erty of invariance under dilation, r — > ar [44[ . The com- 
putation of the first-passage exponent in a generalized 
rf-dimensional cone requires computation of the lowest 
eigenvalue A of the angular portion V 2 of the Laplacian 



V 2 / 



-A/, 



(33) 



subject to the Dirichlet boundary condition, / = 0, on 
the surface of the cone. The choice of the smallest eigen- 
value ensures that / > inside the allowed region. The 
first-passage exponent j3 is related to the lowest eigen- 
value via A = 2/3(2/3 + d - 2), see 

We now invoke a theorem that, in its simplest form, 
states that amongst all domains with the same volume, 
the lowest eigenvalue of the Laplace-Dirichlet operator 
occurs when the domain is a ball. Rayleigh conjectured 
this result for two dimensions [45| , and Faber and Krahn 
eventually proved it [46] . The Rayleigh- Faber-Krahn the- 
orem generalizes to higher dimensions and under mild 
conditions, to Riemannian manifolds [47] . In the context 
of first-passage processes, this theorem implies that of all 
generalized cones of fixed solid angle, the first-passage ex- 
ponent is minimal for the circular cone. Since the regions 
discussed in our study qualify as generalized cones, re- 
sults of the cone approximation constitute rigorous lower 
bounds for the exponents. 

Cones can also provide upper bounds. To establish an 
upper bound, we choose a cone that is inscribed by the 
absorbing boundary. Clearly, the first-passage process is 
faster inside the inscribed cone and therefore, the corre- 
sponding exponent must be an upper bound. 

The cone approximation appears to be asymptotically 
exact for the leadership problem as it predicts the bulk 
of the exponents with excellent accuracy. A first step 
toward proving this exactness is the observation that all 
opening angles in (|26p approach a right angle as the num- 



9 



ber of particles diverges. The other major challenge is 
finding an exact, or at least asymptotically exact, frame- 
work for calculating the family of exponents in the origin- 
crossing problem. 

Finally, we mention that in the absence of a theoretical 
method for obtaining statistics of first-passage exactly, 
Monte Carlo simulations play vital role. We utilized a 



straightforward simulation technique. It is interesting to 
find out how accelerated Monte Carlo simulations fare in 
producing accurate estimates for the exponents [HHHO]. 

This research has been supported by DOE grant DE- 
AC52-06NA25396 and NSF grant CCF-0829541. 
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